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ABSTRACT 

To elucidate the permeation of cosmic ultraviolet (UV) background radia- 
tion into a pregalactic cloud and subsequent ionization, the frequency-dependent 
radiative transfer equation is solved, coupled with the ionization process, for a 
spherical top-hat cloud composed of pure hydrogen. The calculations properly 
involve scattering processes of ionizing photons which originate from radia- 
tive recombination. As a result, it is shown that the self-shielding, although 
it is often disregarded in cosmological hydrodynamic simulations, could start 
to emerge shortly after the maximum expansion stages of density fluctuations. 
Quantitatively, the self-shielding is prominent above a critical number density of 
hydrogen which is given by n crit = 1.4 x lCT 2 cm~ 3 (M/1O 8 M ) -1 ^ 5 if/ 5 for 10 4 K 
gas, where M is the cloud mass and the UV background intensity is assumed 
to be I u = 10~ 21 /2i(^/^L) _1 ergcm -2 s _1 ster -1 Hz -1 with vl being the Lyman 
limit frequency. The weak dependence of n cr i t upon the mass is worth noting. 
The corresponding critical optical depth (r crit ) turns out to be independent of 
either M or I21, which is r cr j t = 2.4 for 10 4 K gas. The present analysis reveals 
that the Stromgren approximation leads to overestimation of the photoioniza- 
tion effects. Also, the self-shielded neutral core is no longer sharply separated 
from surrounding ionized regions; a low but noticeable degree of ionization is 
caused by high energy photons even in the self-shielded core. The present re- 
sults may be substantial on considering the biasing by photoionization against 
low-mass galaxy formation. 
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1. Introduction 



Recently, a recalcitrant problem on galaxy formation has been pointed out that low- 
mass galaxies are overproduced as compared with observations in the context of the hierar- 
chical bottom-up theory of galaxy formation (e.g. White & Frenk 1991; Kauffman, White, 
& Guideroni 1993; Cole et al. 1994). Hence, some process that inhibits the formation of 
low-mass galaxies is required. Photoionization has been considered as one of such mecha- 
nisms (Dekel & Rees 1987; Babul & Rees 1992; Efstathiou 1992; Chiba & Nath 1994; Thoul 
& Weinberg 1996; Quinn, Katz, & Efstathiou 1996). In photoionized media, the cooling 
efficiency is dramatically reduced in the temperature range of 10 4 K < T < 10 5 K. Also, a 
bulk of energy could be carried into a cloud, so that the enhanced thermal pressure could 
suppress the gravitational collapse of a subgalactic cloud with the virial temperature lower 
than several 10 4 K, i.e., M<1O 9 M in gas mass (Umemura & Ikeuchi 1984, 1985; Ikeuchi 
1986; Rees 1986; Bond, Szalay, & Silk 1988; Steinmetz 1995; Thoul & Weinberg 1996). In 
cosmological hydrodynamic simulations, an optically thin medium against ionizing photons 
has been mostly assumed so far (Umemura & Ikeuchi 1984, 1985; Thoul & Weinberg 1996; 
Quinn, Katz, & Efstathiou 1996). In semi-analytic approaches, naive analytic corrections 
for opacity effects have been made based on the optical depth criterion (Efstathiou 1992) 
or the Stromgren approximation (Chiba & Nath 1994). In the case of interstellar clouds, 
it is claimed that the Stromgren approximation could be misleading against the real effect 
of the penetration of diffuse UV (Flannery, Roberge, & Rybicki 1980; Maloney 1993), and 
therefore the radiative transfer equation should be solved properly. However, as far as we 
know, no study hitherto has been made on determining the ionization structure inside a 
pregalactic cloud by solving radiative transfer equation for diffuse UV photons. Hence, the 
effects of photoionization on the evolution of a pregalactic cloud have not been assessed 
satisfactorily. In this paper, we solve the radiative transfer of diffuse UV radiation cou- 
pled with ionization process to elucidate the self-shielding of pregalactic clouds from UV 
background radiation and provide a practical criterion for the self-shielding. 



2. Radiative Transfer with Ionization Process 

We assume for simplicity a spherical top-hat (uniform) density distribution of a cloud, 
and place 100 radial meshes for solving radiation transfer. Also, the cloud is assumed to be 
composed of pure hydrogen so that we could readily compare the results of the frequency- 
dependent radiative transfer with an analytic estimate, although we should keep in mind 
that helium of cosmic abundance could alter the ionization degree maximally by order 10% 
(Osterbrock 1989; Nakamoto et al. 1997). 

The intensity of UV background radiation at high redshifts is inferred from so-called 
proximity effect of Lyman a absorption lines in QSO spectra (Bajtlik, Duncan, & Ostriker 
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1988; Giallongo et al. 1996). The observations require the diffuse UV radiation to be 
at a level of I VLi0 = 10" 21±0 - 5 ergs cm" 2 s" 1 ster" 1 Hz -1 at the hydrogen Lyman edge 
at z = 1.7 — 4.1. In this paper, we assume the specific intensity of UV background as 
I Vj0 = 10 _21 / 2 i(^ '/ 'v L ) n ergs cm" 2 s" 1 ster" 1 Hz" 1 , where we set n — — 1 and vary 7 2 i in 
the range of 0.1 < I21 < 2. As for ionization process, we presume the ionization balance, 
because the ionization or recombination timescale for clouds with density of interest is 
much shorter than the dynamical timescale. The equation of the ionization balance is 



where xhi is the fraction of neutral hydrogen, T 7 is the photoionization rate, T cl is the 
collisional ionization rate, n is the hydrogen number density, and a^(T) is the total recom- 
bination coefficient to all bound levels of hydrogen. T 7 is given by 



with the photoionization cross section being a u = 6.3 x 10 18 (z/^/z/) 3 [cm 2 ], where the local 
UV intensity I v {r) at radius r is determined by solving transfer equation. r c * is given by 
r ci = 1.2 x 10- 8 r 4 1/2 e- 15 - 8 / T4 [cm 3 s- :l ] with T 4 = T/10 4 K. We assume T 4 = 1 except where 
other values are specified. a A {T) is well fit by a A {T) = 2.1 x 1O" 13 T 4 ~ 1/2 0(16/T 4 )[cm 3 s" 1 ], 
where <j>(y) = 0.5(1.7+lny+l/6y) for y > 0.5 or y(-0.3+1.2 \ny)+y 2 (0.5-\ny) for y < 0.5 
(Sherman 1979). 

Photoionization and recombination processes can be respectively regarded as extinc- 
tion and emission with respect to ionizing photons. Thus, the radiative transfer equation 
for ionizing photons is described as 



where Xu is the extinction coefficient (xu = dunxni) and r\ v is the emissivity. If a free 
electron recombines directly to the ground state of hydrogen, the emitted photon has 
enough energy to cause further photoionization. But, when an electron is captured to 
an excited state of hydrogen, the emitted photon does not have enough energy to ionize 
hydrogen, because the kinetic energy of a free electron is of order 0.1 Ryd for ~ 10 4 K 
gas. Thus, the former process is regarded as scatterings which provides the emissivity, 
while the latter process is pure absorption. The effective scattering albedo is given by uo = 
\ola{T') — &b{T)}/c(a{T), which is 0.4 at 10 4 K (e.g. Osterbrock 1989), where ob(T) is the 
recombination coefficient to all excited levels of hydrogen. Hence, we set r] u — for v > i>l, 
and r) VL = hvuaA'n e n p /4:7i5b' for v = u L where 5u = kT/h and n p is the proton number 
density. 

Taking account of frequency-dependence of the emergent UV intensity, it is convenient 
to divide the photoionization rate into two parts: T 7 = T 7 + Tl >UL . As for Lyman limit 



^Xhi + T ct XHI (l - X Hi)n = a A (T)(l - X Hifn, 



(1) 




(2) 



as 



(3) 
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photons, the transfer equation is the integro-differential equation which includes a source 
term by scattering processes of recombination photons. Thus, the equation is numerically 
solved with including iterative procedure. We, without invoking on-the-spot approximation 
which could be misleading for large mean-free-path photons, solve the equation by means 
of an impact parameter method of high accuracy, which allows us to treat diffuse photons 
correctly (Stone, Mihalas, & Norman 1992). We deal with 156 impact parameters for light 
rays. Also, in order to converge the intensity, we employ the Lambda-iteration method 
(e.g. Mihalas & Mihalas 1984). 

In the range of v > since photoionization is regarded as pure absorption, the UV 
intensity is obtained just by l v = I v $e ~ Tv , where I v>0 is the boundary intensity and r v is 
the ionization optical depth. So, TZ> Vh is figured out by 

?i> VL = r dn r ^^r-^ (-) 3 ■ ^ 

L JO Ju L +5u hv V V ) 

Considering that the optical depth is r v = t Ul (v l / vfnxui with the Lyman limit optical 
depth t Vl and that the assumed boundary intensity is I u>0 = l VL ^{y L jv)^ the above equation 
can be analytically integrated using the incomplete gamma function 7 if xhi is given: 

™ _ (*" dc foJvLfl 7(4/3, r VL ) 

With use of this integration, the overall procedure to solve the transfer equation is as 
follows: 

1. Initially, give the cloud mass M and the radius R (therefore the density), and set 
Xhi by assuming optically thin medium. Specify the boundary UV intensity by J21. 

2. Solve the transfer equation at v = v L for given Xhi, and calculate 1^ . Obtain Tl >UL 
analytically by equation (^|). Then, get the total photoionization rate T 7 . 

3. Solve ionization equilibrium using above-obtained T 7 , and thereby renew Xhi- 

4. Continue steps 2 and 3 until xhi converges at a level of relative error of 10 -6 . (Typ- 
ically 100 iterations are performed.) 

Note that the analytical integration @ with respect to frequencies enables us to 
reduce the computational cost dramatically on solving such a frequency-dependent transfer 
equation including scatterings. The validity of this method is confirmed by exactly solving 
the transfer equation with using a number of meshes for frequencies. 
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3. 



Numerical Results 



We consider the cloud mass range of M — 10 5 ~ 9 M Q and vary the radius in the range 
of R = 0.1 — 25kpc. Fig.l shows the growth of self-shielded regions when a cloud of 1O 8 M 
is contracting, embedded in the UV background of I21 — 1. The self-shielding is prominent 
when i?<4kpc. It is noted that the distributions of the HI fraction are a gradual function 
of radii even in the self-shielded stage, where we cannot recognize a clear boundary between 
the neutral core and the ionized envelop, and also a low but noticeable ionization is left 
in the self-shielded regions. Such distributions seem to be realized by the high frequency 
photons far above the Lyman edge which could permeate into deeper regions due to the 
smaller cross section for ionization. To ensure this conjecture, we solve the transfer solely 
for the Lyman edge photons to obtain I VL (r) and tentatively set the form of UV radiation 
spectrum to be I v (r) = I UL {r){vi/ v) at any radius (which implies that photons of higher 
frequencies are absorbed with the same ionization cross section as that at the Lyman 
edge). We can see an outstanding difference between two cases as shown in Figure 1. In 
the tentative case, we can see a steep inward increase of neutral fraction and a very sharp 
transition from ionized regions to neutral regions. 

Figure 2 shows the neutral hydrogen fraction at the center as a function of cloud size. 
We see that the central xhi varies abruptly at a certain critical size. Here, we define the 
critical radius (R cr it) at which the HI fraction at the center drops just below 0.1. In Figure 
3, we plot the critical radii obtained from all the numerical results as a function of UV 
background intensity. Different symbols represent different mass of gas clouds. All the 
results can be remarkably well fitted by a simple formula which is a function of the cloud 
mass and the UV background intensity; 



It is worth noting that the critical density is quite weakly dependent upon M. If the cloud is 
highly ionized and optically thin, the neutral fraction should be Xhi,o — O.lSn/^ 1 with ne- 
glecting collisional ionization. Then, we define the critical optical depth r cr j t at the Lyman 
edge as a measure such that the self-shielding becomes effective: r cr n = n cr i t XHifl a 'v L Rcrit — 
2.4, which turns out to be independent of not only M but also 7 2 i- Hence, such a simple 
criterion of r = 1 as adopted by Efstathiou (1992) is found heuristically to be practical in 
order to assess the self-shielding for ~ 10 4 K clouds. 

For clouds of different temperature, we have found that R cr u is scaled by as (T) 1 / 5 , 
and therefore n crit oc chb(T)~ 3 / 5 . Resultingly, r crit is scaled by a j4 (T)/Q; B (T)[= (1 — a;) -1 ]. 





Equivalently, the corresponding critical number density of the cloud is 
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Both Qa{T) and a_e(T) are decreasing functions of temperature, but uj increases with 
temperature, so that r cr u is larger for higher temperature. For instance, T cr u = 2.7 for 
3 x 10 4 K. In an extreme case of infinite temperature (although unrealistic), the complete 
scattering {u = 1) leads to r crit = oo. In other words, clouds with any optical depth can 
be ionized due to photon diffusion. 



4. Comparison with Analytic Estimates 



Here, we try to analytically estimate the critical radius based upon the Stromgren ap- 
proximation. (A similar estimate is found in Chiba & Nath 1994.) By equating the number 
per unit time of ionizing photons which enter from the surface to the number per unit time 
of photons which are absorbed in the cloud, we have Rhi = [R 3 — 3'irR 2 I 1/Lt o/hn 2 (l — Xhi) 2 
as(T)] 1 / 3 , where the ionized regions are assumed to be sharply separated from the neutral 
core of radius Rhi- Then, the critical radius can be estimated by setting Rhi = 0: Resul- 
tantly we find R CT u = 3.5kpc (M/1O 8 M ) 2//5 I^i^ , and the corresponding critical optical 
depth is r crit = 5.3. Hence, the dependence of @ upon the cloud mass and the UV inten- 
sity can be fundamentally understood by this argument. But, from a quantitative point 
of view, this approximation obviously leads to overestimation of photoionization effects as 
recognized by the critical optical depth. The overestimation comes from the assumption 
that all photons which enter into the cloud always cause ionization. In fact, some photons 
which especially have low incident angles do escape from the gas cloud without causing 
ionization. The diffusion process of ionizing photons tends to enhance this effect. Fur- 
thermore, a sharply edged neutral core, which is the basic assumption in the Stromgren 
approximation, is no longer realistic as shown above. 



5. Discussion 

The present results seem of a great significance on considering the biasing by photoion- 
ization against the formation of low-mass galaxies. It is shown in previous analyses that if a 
cloud is assumed to be optically thin, the photoionization suppresses the collapse of clouds 
with M < 1O 9 M (Umemura & Ikeuchi 1984) or circular velocities smaller than 30 km s _1 
(Thoul & Weinberg 1996). In the present analysis, the permeation of UV radiation is 
characterized by a different criterion, and the critical density has turned out to be almost 
independent of the mass. Hence, the evolution of subgalactic clouds would not be deter- 
mined solely by the cloud mass or the circular velocity. The maximum expansion radius of a 
top-hat density fluctuation is given by R ma x = 10.7 kpc(M/10 8 M o ) 1//3 [10/(l + z max )}h^^ 3 
in an Einstein-de Sitter universe, where z max is the maximum expansion epoch, /i 50 is the 
present Hubble constant in units of 50km s _1 Mpc _1 , and the baryon density parameter is 
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assumed to be 0.05. Comparing R max with (|6|) and taking into account a possibility that the 
UV intensity might be significantly lower at z > 5, we can speculate that the self-shielding 
could be quite effective shortly after the maximum expansion stages. Thus, in order to 
assess the effects of photoionization properly, we should consider the frequency-dependent 
radiative transfer of diffuse UV photons. 

The present results are also relevant to the formation of primordial hydrogen molecules 
which provide key physics for the formation of the first generation objects (Tegmak et al. 
1997). UV radiation will naively suppress the formation of hydrogen molecules and thereby 
cooling (Haiman, Rees, & Loeb 1997). However, when the reionized gas is self-shielded in 
the course of evolution, H~ or H-2~ ions could form efficiently due to the residual ionization 
which may lead to the effective production of hydrogen molecules (e.g. Kang & Shapiro 
1992). Since the formation of primordial hydrogen molecules is a non-equilibrium process, 
the permeation of even a small portion of UV background photons may play an important 
role for subsequent molecule formation, and thereby the evolution of neutral core. 

We are grateful to T. Nakamoto, and H. Susa for helpful discussion. This work was 
carried out at the Center for Computational Physics of University of Tsukuba. This work 
was supported in part by the Grants-in Aid of the Ministry of Education, Science, and 
Culture, 09874055. 

REFERENCES 

Babul, A., & Rees, M. J. 1992, MNRAS, 255, 346 
Bajtlik, S. Duncan, R. C. & Ostriker, J. P. 1988, ApJ, 327, 570 
Bond, J. R., Szalay, A. S., & Silk, J. 1988, ApJ, 324, 627 
Chiba, M., & Nath, B. 1994, ApJ, 436, 618 

Cole, S., Aragon-Salamanca, A., Frenk, C. S., Navarro, J. F., & Zepf, S. E. 1994, MNRAS, 
271, 744 

Dekel, A., & Rees, M. J. 1987, Nature, 326, 455 
Efstathiou, G. 1992, MNRAS, 256, 43 

Flannery, B. P., Roberge, W., & Rybicki, G. B. 1980, ApJ, 236, 598 

Giallongo, E., Cristiani, S., D'Odorico, S., Fontana, A., & Savaglio, S. 1996, ApJ, 466, 46 

Haiman, Z., Rees, M. J., & Loeb, A., 1997, ApJ, 476, 458 

Ikeuchi, S. 1986, Ap&SS, 118, 509 

Kang, H., & Shapiro, P. R. 1992, ApJ, 386, 432 



-8 - 



Kaufman, G., White, S. D., & Guideroni, B. 1993, MNRAS, 264, 201 
Maloney, P. 1993, ApJ, 414, 41 

Mihalas, D., & Mihalas, B. W. 1984, Foundation of Radiation Hydrodynamics (New York: 
Oxford Univ. Press), 366-369 

Nakamoto, T., Susa, H., & Umemura, M. 1997, in Proc. International Symposium on 
Supercomputing, New Horizon of Computational Science, in press 

Osterbrock, D. E. 1989, in Astrophysics of Gaseous Nebulae and Active Galactic Nuclei 
(University Science Books) 

Quinn, T., Katz, N., & Efstathiou, G. 1996, MNRAS, 278, L49 

Rees, M. J. 1986, MNRAS, 218, 25p 

Sherman, R.D. 1979, ApJ, 232, 1 

Steinmetz, M. 1995, in Proc. 17th Texas Symp., Relativistic Astrophysics, ed. H. Bohringer, 
G. E. Morfill, & J. E. Triimper (Ann. NY Acad. Sci., Vol. 759), 628 

Stone, J. M., Mihalas, D., & Norman, M. 1992, ApJS, 80, 819 

Tegmak, M., Silk, J., Rees, M. J., Blanchard, A., Abel, T., & Palla, F. 1997, ApJ, 474, 1 
Thoul, A. A., & Weinberg, D. H. 1996, ApJ, 465, 608 
Umemura, M., & Ikeuchi, S. 1984, Prog. Theor. Phys., 72, 47 
Umemura, M., & Ikeuchi, S. 1985, ApJ, 299, 583 
White, S. D. M., & Frenk, C. S. 1991, ApJ, 379, 52 



This preprint was prepared with the AAS IATgX macros v4.0. 



-9 - 



Fig. 1. — The ionization structure is shown by the neutral hydrogen fraction xhi as a 
function of radii in units of the cloud radius (thick lines). Here, M = 1O 8 M and J 2 i = 1 
are assumed. It is seen that the low ionization regions rapidly grows due to self-shielding 
when the radius is smaller than 4.4kpc. It is a characteristic feature that the boundary 
between the neutral core and the ionized envelop is not clear in contrast with the Stromgren 
sphere. This is because a hard power-law type (n = —1) of UV background radiation 
causes appreciable ionization even in deeper regions by high frequency photons far above 
the Lyman limit. Consequently, a small ionization fraction is left in the neutral core. 
Comparatively, if an artificial softer spectrum of UV is assumed, a very narrow transition 
region appears and the ionization in the neutral core is negligible (thin lines). 

Fig. 2. — The central neutral hydrogen fraction is shown against the cloud radius for a 
wide variety of cloud mass. /21 = 1 is assumed here. It is shown that the self-shielding is 
prominent below a critical radius (which is defined by the condition that xhi = 0.1 in the 
present paper). 

Fig. 3. — The critical radii are plotted as a function of UV background intensity. Different 
symbols represent different mass of gas clouds. The straight lines are the fitting formula 
(see the text). 
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